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Abstract 

Several aspects of next-to-leading (NLO) order corrections to see-saw formulas are 
discussed and phenomenologically relevant situations are identified. We generalize 
the formalism to calculate the NLO terms developed for the type I see-saw to variants 
like the inverse, double or linear see-saw, i.e., to cases in which more than two mass 
scales are present. In the standard type I case with very heavy fermion singlets the 
sub-leading terms are negligible. However, effects in the percent regime are possible 
when sub-matrices of the complete neutral fermion mass matrix obey a moderate 
hierarchy, e.g. weak scale and TeV scale. Examples are cancellations of large terms 
leading to small neutrino masses, or inverse see-saw scenarios. We furthermore iden- 
tify situations in which no NLO corrections to certain observables arise, namely for 
fi~r symmetry and cases with a vanishing neutrino mass. Finally, we emphasize that 
the unavoidable unitarity violation in see-saw scenarios with extra fermions can be 
calculated with the formalism in a straightforward manner. 
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1 Introduction 



Neutrino masses are small. This simple fact is usually attributed to the presence of a see- 
saw mechanism. In its simplest and most often studied manifestation, the type I see-saw 
PQ, the low energy Majorana neutrino mass matrix is 

m v = -m T D M^ 1 m D , 

where is a Dirac and Mr a Majorana mass matrix. The above relation is formally 
obtained in the limit of "Mr ^> mp", which means that the eigenvalues of Mr are much 
larger than the entries of nip. While one usually gives the expression for m u with a equality 
sign "=", it should strictly speaking be an approximative equal sign as there are 

higher order, or next-to-leading order (NLO), terms which correct it. In the standard type 
I see-saw with heavy fermion singlets the corrections are usually negligible. In the present 
paper we will study the NLO corrections in detail and identify phenomenologically relevant 
applications. A formalism to give those terms at arbitrary order has been developed by 
Grimus and Lavoura in Ref . [2] . We generalize that formalism^ to several of the see-saw 
variants which are discussed in the literature, such as the double [I], inverse [5], linear [6] 
or singular [7] see-saw. We note that in scenarios in which the see-saw scale is lowered to 
TeV scale, NLO corrections in the percent regime are possible, which in the light of up- 
coming neutrino precision experiments is surely not negligible^. This occurs for instance in 
scenarios in which cancellations of large terms lead to small neutrino mass, or in inverse see- 
saw frameworks. These cases have in common that sub-matrices of the complete neutral 
fermion mass matrix obey a moderate hierarchy. 

One common aspect of all see-saw mechanisms with additional fermions is the violation 
of unitarity of the 3x3 mixing matrix which describes the mixing of the three active 
neutrinos. The formalism to calculate NLO see-saw corrections is applicable to calculate 
the magnitude and structure of those. This allows to obtain formulae for unitarity violation 
in a simple and straightforward manner for the see-saw variants. 

We furthermore identify situations in which no corrections to certain parameters arise. We 
will show that in case of /x-r symmetry there are no NLO corrections to U e ^, = and 
$23 = vr/4. Another finding is that if one neutrino is massless, then the mixing matrix 
elements associated with this massless state receive no corrections. The massless neutrino 
does not mix with the heavy ones. 

The paper is build up as follows: in Section [2] the formalism to calculate higher order 
corrections to the type I see-saw term is reviewed. The connection to the inherent unitarity 
violation is noted. In Section [3] we discuss examples on the application of the NLO terms, 
and identify cases which are stable. In Section [4] we show how the formalism can be applied 
to see-saw variants, before we conclude in Section [5j 

*A similar Ansatz for the unitary matrix diagonalizing the full neutral fermion mass matrix has been 
proposed in Ref. [3]. The same results could be obtained with this approach. 

^We will discuss here only the corrections to the low energy mass matrix m u and leave consequences 
for the heavy singlets for further study. 
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2 NLO Terms to the Type I See-Saw Mechanism 



In what follows we will review the derivation of the NLO terms to the type I see-saw for- 
mula. The reader familiar with it can continue in Section l2~2l where some of its applications 
are studied. 



2.1 Derivation of NLO Terms 

The conventional type I see-saw mechanism pQ contains after electroweak symmetry break- 
ing two mass terms in the Lagrangian: 

£ = N^m D u L + ^M R N C R + h.c. (2.1) 

Here Ul are left-handed neutrinos, Nr are right-handed singlets, and mjj (Mr) is the 
Dirac (Majorana) mass matrix. The effective mass matrix at low energy is conventionally 
obtained by integrating out the heavy states Nr, and the result is 



rriv 



-m T D M^ 1 m D . (2.2) 



The approximative nature of this expression is noteworthy. Formally, the effective mass 
matrix is obtained by diagonalizing the total mass matrix in the basis (z/£, Nr): 

We note here for later use that the inverse of a matrix of this texture is 

ka-1 _ ( - m ~D MRiml)' 1 m D x \ 

M "V (O' 1 o )■ (2 - 4) 

We will assume in what follows that the involved matrices are invertible square matrices, 
unless otherwise noted. However, in cases with non- invertible (Sections 13.21 and 13 .3 j) 
we will be able to exactly solve the problem without the need of an expansion. If the 
eigenvalues of Mr are all much heavier than the entries of m^, then this situation will be 
described as U M R ^> m D " . Block diagonalization of M. is now possible, and the block- 
diagonal matrices are approximately given by 

-m T D Mn l m D and Mr. (2.5) 

In Ref. [2] a formalism to evaluate the corrections to these expressions to arbitrary order 
has been developed. Let us shortly summarize the derivation of that result. A unitary 
transformation diagonalizes Ai according to 
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and transforms the states (i>l,N r ) to the mass states where the subscript "1" 
denotes light and "h" denotes heavy: 

"UK:).- (2 - 7) 

The matrix 14 can be written as [2] 



U ~\ -fit vT^Wb ) ' U - \ fit VT^Wb ) ' (2 ' 8) 

where B is a complex 3x3 matrix (in general it has the dimension of mo), and the square 
root is to be understood as 

\fl - BB^ = 1 - - RB+ - - BB ] BB' - . . . - ^Olt^. (BB ] ) n (2.9) 
V 2 8 n\T{-\) y ' K ' 

With this Ansatz the matrix U is unitary order by order in BB^ . The analogy of the form 
of U with a real two-by-two mixing matrix is obvious. We can insert U in Eq. (I2.6P and 
the result for the three independent entries of the r.h.s. is 

y/l - B*B T m T D ^1-BW - B* m D B - B* M R y/l - B^B = , 

-B* m D yjl - BB^ - y/l - B*B T m T D B^ + B* M R B^ = m u , (2.10) 

Vl - B T B*m D B + B T m T D y/l - B^B + \J \ - B T B* M R y/l - B^B = M R . 

Now the see-saw approximation enters the game, by assuming that B can be written as a 
power series in terms of 1/M R , i.e., in terms of the eigenvalues of M R , which are assumed 
to be much heavier than the entries of mo. Hence, B = B\ + B 2 + . . ., where -B, is of order 
(1/M R ) 1 . The square root then reads 

\A - BTB B[B X - l - [b[B 2 + B\B^ (2.11) 

A recursive solution of Eq. ( I2.10p is now possible. At leading order the solution of Eq. ( I2.10p 
is given by 

B\ = m T D M- 1 . (2.12) 

The next order term B<i in the expansion is, in the limit of a vanishing triplet contribution, 
zero. This is true for all Bi, where % is even [2]. We obtain for the 3rd and 5th order terms 

B; M r = - l -BlB T v m r D - 1 - m T D B\ B 1 - B{ m D B 1 + ~ B{ M R B\ B x 

= -\m T D M- 1 (M* y 1 m* D m T D - m T D M" 1 m D m^ D (M* )~ x , (2.13) 

b; m r = - i (bi bI + b; b t x + ~ b\ bJ b\ b^J m T D - b\ m D b 3 

-B*m D B 1 + ^B*M r B\B 1 . 
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Inserting B\ and B% in the 11- and 22-entries of Eq. (12. 6 j) yields 



m u = — m D M R mo + - M r m D (M 



(M 



D " b D 



M R X m D 



M R =M r + - 



r*\— 1 * T 

m D m D 



m D ml (M*)- 1 + (M* ; 

One is lead to define the symmetric matrix 

X = A + A T , where A = m D m ] D (M^)" 1 . 
The order of magnitude of X is m 2 D /Mn and one can simplify the relations to 



(2.14) 



m. 



M 



R 



- m T D M^ 1 m D + -m T D M R l X M' 1 m D 
M R + l -X. 



For completeness, we also give the NNLO terms to m u and Mr, which are 



-NNLO 



m 



--r:m T D M^ 



(2.15a) 
(2.15b) 



\am^a + \a t m^a t + \a t m^a+ 1 -{m* r )-^a*a t 



+-AA^ (Mr)" 1 + AA* (Mr)" 1 + (M*rY 1 A ] A 
2 



T 



M~ x m D 



M^NLO 



A A* {M* R )- 1 + (Mr .y 1 A^A T + ^A M^ 1 A + ^A T M R l A T 



The zeroth order terms of the light and heavy mass matrices are m 2 D /MR and Mr, respec- 
tively. The relative NLO corrections are of order X/Mr = rn? D jM\ for both. The absolute 
correction is of order mjj/M R for the light neutrinos and m 2 D /MR for the heavy neutrinos. 
Note that this NLO correction vanishes when A is antisymmetric. The absolute order of 
magnitude of the NNLO terms is m%/M R for the light neutrinos and m A D /M R for the heavy 
ones. In general, the N n+1 LO term of the heavy neutrinos has the same absolute order of 
magnitude than the N"LO term of the light ones. 



A comment to be made here is that the same expressions for the corrections are obtained 
in type III see-saw scenarios [IT] , for which Mr is the mass term of the neutral component 
of a weak fermion triplet. 

In this work we will mostly ignore the possibility of the presence of a Higgs triplet, which 
would fill the upper left entry of A4 in Eq. (12.31) with a term [10]. In this case, the first 
order correction to m v is [2] 

m v = m L -m T D M^m D + )^m T D Ml x XM R l m D - l -(C + C T ) , (2.16) 

where C = mF D M R l (M* R )~ l m* D m L . 
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2.2 Some possible Applications 

We will continue with a few examples on the possible consequences of the NLO terms. 
The typical order of magnitude of the terms is vtld ~ 10 2 GeV and Mr ~ 10 14 GeV, for 
which m v ~ 0.1 eV. In this case, of course, the NLO terms are negligible (the same is true 
for the unitarity violation, see below). An exception is when the Majorana singlets are 
put to TeV scale, which is largely motivated by current collider opportunities. The mixing 
with the singlets is naively of order mo /Mr, and hence the requirement of m v = 0.1 eV 
would lead to small mo and thus small mixing. However, it is possible that small neutrino 
masses are an effect of cancellation of large terms, with Mr ~ TeV and mo ~ v , in which 
case the ratio of leading order and NLO terms is m 2 D fM\ ~ 10~ 2 , a percent effect! In 
the light of future precision experiments, this is a correction one surely should take into 
account. 

Let us give an illustrative example on this. We enter now the basis in which Mr is real 
and diagonal. For a two neutrino case, and a diagonal right-handed neutrino mass matrix 
Mr, and if the Dirac mass matrix is written in its most general form 

m D = ( a c b d ) , (2.17) 

the neutrino mass matrix is at leading order: 

1 / a 2 ab \ 1 / c 2 cd \ , n m\ 

m " = w\ ■ « ) + wA ■ * )■ i2M > 

One might imagine that M\^ lie around TeV, ttld lies around v, and that at leading order 
the neutrino mass matrix vanishes, m v = 0. In order for m u to vanish^, the requirements 

m D =( X aX ) and ^ = (2.19) 
\V ay J y 2 M 2 V ' 

must hold simultaneously [12] . A small correction in one or both of those two conditions 
will generate small but non-zero neutrino mass. For instance, violating the second condition 

as x = iy a ^ (1 + e), one finds 



ml = -mlM R l m D = ^-e(2 + e)(^ 1 °Y (2.20) 

which has eigenvalues mi = and m 2 = (1 + a 2 ) y 2 /M 2 (2 + e) e, obtained by diagonalizing 
with the mixing matrix 



-a 1 



VI + a 2 V 1 a 



tThis will hold true to all orders, see below. 
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Note that large mixing would imply a = 0(1). For, a = 1, y ~ 100 GeV and M 2 — TeV, 
it follows that e ~ 10~ 12 to give m 2 = 0.05 eV. This illustrates the enormous tuning which 
has to present in order for this mechanism to work. 

What are the corrections to m"? We have shown above that the NLO term to m" is 

jr, (M^ 1 )*. Evaluating this with 



\m T D M^ 1 X M^ 1 m D , where X = A + A T and A 



rriD ra 



our example gives 



m. 



(2 



e e 



M l M% 



(M 1 + (l + e) 2 M 2 ) y 4 (l 



a 

,:2 



(2.21) 



Comparing the zeroth and first order term, we havem° = 0(ey 2 /M) and ml = 0(ey 4 /M 3 ), 
which for y ~ 100 GeV and M ~ TeV results in a NLO term being suppressed only at 
the percent level. The mixing matrix stays the same in this example, and the non-zero 
eigenvalue of m° v + ml is different from the non-zero eigenvalue of m° by order ey A /M 3 . 
It should be clear to realize that the moderate hierarchy between the submatrices mo and 
Mr of Ai leads here to sizable effects. 



One may wonder what happens when m v 



m D 




at leading order, for instance if 


and Mr 




(2.22) 



However, one can show [HI [2] that the vanishing of m v will remain true to all order^j] (the 
rank of M. is three). We will discuss the realistic case of one vanishing eigenvalue of m v 
in Section l3~3l 



Another aspect is that zeros entries could be filled by NLO terms. Consider the case 



ci 
me =| 62 c 2 I and Mr 

a 3 b 3 




(2.23) 



The resulting low energy mass matrix at zeroth order is 

m T D M^ 1 m D 



' U AI 2 M 2 * 

9 b2 &3 b 3 c 2 
* M 2 M2 



V 



(2.24) 



Ah 



h J 



and forbids neutrino-less double beta decay because the 11-entry of m v vanishes. The NLO 
term fills this entry with a contribution 



/ % a\ b 2 b 3 
(^)n =- 



Ml ■ 



(2.25) 



§ A proof for the case when a triplet is present can be found in [5J. 
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However, this term is suppressed with respect to the non-zero entries in mF D M^ 1 rrijj by 
many orders of magnitude, and can safely be neglected. It is typically even smaller than 
terms of order U^mf/q 2 (q 2 ~ 0.1 GeV 2 the typical momentum exchange in double beta 
decay), which are in general non-zero when U 2 ^ rrii = 0. An exception could be when such 
zero textures are generated in scenarios in which small m v is generated by cancellations of 
large terms, see above. 



2.3 Connection to Unitarity Violation 

It is well-known that an intrinsic unitarity violation is present in those see-saw scenarios 
which contain extra fermions (i.e., not in a pure type II see-saw). This can be shown easily 
by diagonalizing Eq. (12. 6 p in a slightly different way, namely via 

w=(r v)' (2 ' 26) 

where N, S, T, V are 3x3 mixing matrices which are in general non-unitary. By evaluating 

K r MM={ W \ > M lg)- (2.27) 

and by assuming that S, T are of order mc/Mfl, one obtains from the 12-entry of U T M. V U 
that at leading order T T ~ — iV T m^M^ . Inserting this in the 11-entry gives 

m diag _ _ N T m T M -l ^ N ( 2 2 g) 

and V T M R V ~ M^ iag . The PMNS matrix N is therefore non-unitary, because NN^ = 
1 — SS^ and N'N = t — T^T. Phenomenologically, the non-unitarity can be described by 
writing 

N = (l + ri)U , (2.29) 

where Uq is unitary and r) hermitian. The latter matrix contains three phases, but current 
constraints exist only for its absolute values [T3j : 



4.0 x 10~ 3 1.2 x 10~ 5 3.2 x 10" 3 

\rj\ < | • 1.6 x 10~ 3 2.1 x 10" 3 | . (2.30) 

5.3 x 10~ 3 

By comparing Eq. (I2.29p with (12.261) we can identify r) ~ —-^SS*. By inserting in 
N*N = 1 — T^T the above relation for T and the phenomenological description for N 
from Eq. (I2.29p . one finds 

V ~ - 1 - rrJ D (M^ 1 )* M^ 1 m D = ~ Btfl . (2.31) 



S 



Hence, we can read off the amount of unitarity violation from the first order expression 
of B, which is given as B 1 = m D (M^)^ 1 in Eq. ( 12.12ft If a triplet term mj, is present, 
the same calculation can be performed and the result from Eq. ( 12.311) stays the same as 
long as mi <C Mr. A triplet term does therefore not induce unitarity violation. As we 
will see in Section HI we can calculate the NLO terms for see-saw variants like the inverse 
or double see-saw in the same way as we have done above for the type I see-saw. This 
makes it possible to simply write down the magnitude and structure of unitarity violation 
for those scenarios. 

3 Special Cases in Type I See-Saw 

There are certain cases in which the specific (flavor) structure of the mass matrices leaves 
imprints on the higher order see-saw corrections and the unitarity violating parameters. 
We will discuss some examples, starting first with /i-r symmetry and then move on to 
scaling (7713 = #13 = in the inverted hierarchy), which we will generalize to scenarios 
containing a vanishing neutrino mass with arbitrary mixing and mass ordering. 

3.1 fi—r Symmetric See-Saw 

Consider fj,-r symmetric!] mo and Mr, i.e., [H] 



As a result of such a structure the leading term of the low energy mass matrix is 



where A,B,D,E are functions of the parameters in and Mr. The above matrix 
predicts to the eigenvalue D — E the eigenvector (0, —1, 1) T . Hence, if D — E corresponds to 
the largest (smallest) mass, #13 = and 6*23 = it/4 in the normal (inverted) mass ordering 

^Actually, the above limits on r\ often assume that the physics leading to unitarity violation is inac- 
cessible at low energy. For instance, the amplitudes of lepton flavor violating charged lepton transition 
li — > £j receive contributions from fermion singlets in the form of Sfc^jfc^ « 9{ x k)> where Xk = m^/m^ 
in the loop function g(x), and rrik is a singlet fermion which mixes with the SM particles. For singlet 
masses above a few 100 GeV the limits are basically equivalent to (|2.30p . 

"Actually, this is a generalized form of the usually considered /x-t symmetry, which denotes the invari- 
ance under exchange of flavor indices fi and r in m„. A more correct name would be 2-3 symmetry, but 
we stick to the name /i— r symmetry. 






(3.2) 
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is predicted. Interestingly, the unitarity violating parameter 7] is also /x-r symmetric, 



— m 



D 




(3.3) 



where the new parameters x, y, z, w are functions of the entries in and Mr in Eq. (13. ip . 
This implies in particular that i] eT is predicted to be extremely small and below values 
which can be probed in future neutrino oscillation facilities |15j . 

The eigenvalue D — E = (e — f) 2 /(w — z) of the zeroth order matrix is in its exact form 
(that is, by diagonalizing M. instead of mj) given as 



(z-w)- ^(z-w) 2 + A(e- /) 2 ) = M - VM 2 + A 
The exact eigenvector to this eigenvalue can be written as 



6 2 



77T 



(3.4) 



/ 


u e3 
























Un 2 3 






Un 3 3 


J 



( 



N 





-1 
1 




\ 



_ J_ ( M - VM 2 + 4m 2 ) 
V ^ [M — a/ M 2 + 4 m 2 ) y 



(3.5) 



where we included its normalization in the factor N. We therefore showed that \J e3 = 
and maximal mixing in the sense {U^l = \U t3 \ is not modified by higher order corrections 
in \x-r symmetric see-saw scenarios. Interestingly, the state with mass m 3 does not mix 
with one heavy neutrino and mixes in equal amounts with the other two. 



3.2 Scaling 

The next example deals with "scaling" [16] (see also [17] ) , for which the Dirac mass matrix 
has the following texture: 

/ a x h bi/c \ 
m D = a 2 b 2 b 2 /c . (3.6) 
\ a 3 b 3 b 3 /c J 

The third column is proportional to the second one, the relevant factor being the "scaling 
constant" c. Interestingly, independent on the form of Mr (other than being non-singular) 
the low energy mass matrix has the form pi 



A 




m 1 D M-'m D =\ ■ D D/c . (3.7) 



Such a low energy mass matrix has been derived for instance in explicit flavor symmetry 
models based on in Ref. ([IS]). The prediction of this particular texture is that the 
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eigenvector to the zero eigenvalue (note that the rank of mF D M^ 1 tud is 2) is (0, — 1/c, 1) T , 
and hence scaling predicts an inverted hierarchy, with U e % = and tan 2 823 = l/|c| 2 [T6] . 
It is easy to see that with mp given in Eq. ( 13. 61) the full 6x6 mass matrix Ai has rank 5 
and the eigenvector corresponding to the zero eigenvalue is 



/ 


u e3 


\ 










u r3 












Un 2 3 




\ 


Un 3 3 


J 



l + IVc| 



/ \ 

-1/c 

1 




V / 



(3. 



1 2 • 
c m 



Therefore, there are no corrections to the predictions U e 3 = and tan 2 623 = 1, 
see-saw scenarios obeying scaling. The massless neutrino does not mix with the heavy 
ones. The unitarity violation obeys the relations 







VfMfJ, 






Ver 




Vfir 




Vtt 



cot 9 



23 



(3.9) 



and again the implied value of r\ er is very small. 



3.3 Vanishing Eigenvalue 



The specific example discussed in the last subsection had a vanishing neutrino mass in 
mP u = —rn^M^mf), and higher order corrections did not induce a non-zero mass (this 
is actually trivial, since the rank of Ai is five), nor did they modify the mixing matrix 
elements of the vanishing eigenvalue. Is this true in general? In what follows we will show 
that this is indeed the case. 

If m° is to have rank 2, it follows that there is an eigenvector l^) to it such that = 0. 

If Mr is non-singular (we treat this case of det Mr = later) then this means that mo has 
rank 2. Hence, there is an eigenvector |$) to mo with the property |$) = 0. With the 
definition of mJJ it follows that m" |$) =0. Since m" can not have two zero eigenvalues, 
|$) must be proportional to and hence mo \ ip) — 0, or to be more specific: 



orra/) 




(3.10) 



depending on whether the normal or inverted ordering is present. Note again that this is 
independent on the form of Mr. Taking the normal ordering as an example, Eq. ( 13 . 1 j) 
implies that the Dirac mass matrix takes the form 







h 


m D = 


a 2 






\ 0,3 


h 



Uel ai+Ufj,i bi \ 
Url 

U el a2+t/ M i b 2 

Url 

Uel Q3+(V)J1 H 

Url J 



(3.11) 
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or similar relations in the first or second column of mp. In case of an inverted hierarchy, 
U a i has to be replaced with U a s, and if in this case U e z = 0, we have the scaling scenario 
described above in Section 13.21 Inserting Eq. (13.111) in the 6x6 matrix Ai reveals that it 
has rank 5, and the exact eigenvector to the zero mass state is simply 



u ei 


\ 























J 



with i — 1 (3) for the normal (inverted) ordering. Thus, if there is a vanishing eigen- 
value rrii of m°, then there are no corrections arising to its mixing parameters U a i, where 
a = e, /i, r. In addition, this neutrino does not mix with heavy ones. 



Is this conclusion valid in both ways? Let us assume that the 6x6 mass matrix Ai has a 
vanishing eigenvalue, i.e., 

Ma = 0, (3.13) 

where a = (ai, a 2 , a%, a^, as, ao) T is the eigenvector of the zero eigenvalue. Solving this 
equation for, say, the third column of mo and 04,5,6 gives nothing but Eqs. A3. 1 lj) and 
( 13.121) . Hence, we have shown that the eigenvector to a vanishing neutrino mass receives 
no corrections from higher order see-saw terms. 



3.4 Almost vanishing m v 

Let us return to the scenarios in which small m„ is generated by a small perturbation to 
m v = 0. In a 3 family framework the condition for vanishing m u is for diagonal Mr that 

(x ax bx 
V ay by 
z az bz 

The generalization to non-diagonal Mr has recently been discussed in Ref. [20], and the 
possibility of percent effects of the NLO terms has been discussed in Section 12.21 The 
structure of the implied unitarity violation has been analyzed recently in Ref. [21]. With 
Eq. (I3.14p and the definition of 77 in terms of B\ we have 

--U^ + 3(::iH 6 )' (315) 

As \r\ ep \ is known to be very small, a can essentially be set to zero and the flavor structure 
of t] becomes very simple |21j . 



x^ z^ 
and M- + M- + M, 



0. 



(3.14) 
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4 NLO Terms to See-Saw Variants 



There are popular variants of the type I see-saw, in which additional singlets S are added 
to the theory, and the (z/£, Nr) basis is extended to a (i/£, Nr, S) basis: 

/ m L m T D m T DS \ 
M = I m D M fi . (4.1) 

\ m DS rriRs M s ) 

The diagonal entries are complex symmetric, while the off-diagonal elements are arbitrary 
complex matrices. As mentioned above, we will not consider the presence of a triplet term 
rriL here. The frequently discussed variants of the type I see-saw are obtained from this 
equation by setting some terms to zero and assuming a hierarchy in the eigenvalues of the 
surviving terms. We will discuss in the following these variants and apply the formalism 
discussed in Section [2]to analyze the NLO terms and the order of magnitude of the unitarity 
violation. The results of this Section are summarized in Table [TJ 



4.1 Double See-Saw 

In the double see-saw scenario we have [3]: 

m T D 
M= | m D m T RS 
m RS M s 

with the conditions m^m^s <C Ms and mo < m RS /Ms- To block-diagonalize Ai, we 
define 



IVJLd 

and write 



m ° ^ and M* := ( ° ) , (4.3) 



a - \m RS M t 



is 



m ={m d m r )- < 4 - 4 ) 

The eigenvalues of the symmetric block M.r are of order Ms and m 2 RS /Ms which are, 
because of the above mentioned conditions, much bigger than the entries in Mp. Thus, if 
we compare Eq. (14. 4 p with Eq. (12. 3p . we recognize that the double see-saw formulas and 
their corrections simply follow from the type I equations which we presented in Section [2j 
We only have to perform the substitutions 

m D ->■ M D , Mr ->■ Mr . (4.5) 

Note that M.r has the same structure as the type I see-saw matrix (12. 3p whose inverse 
is given in Eq. ( 12. 4p . and the inverse of M.r can therefore simply be read off from that 
expression. For illustration, let us determine the usual double see-saw formula and its first 
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order correction. The relevant relations are given in Eq. (I2.15p . Applying (14. 5p . equation 
(I2.14p translates into 

X = A + A T , where A = M D M 1 D (M^)" 1 , (4.6) 
and we obtain for the mass of the lightest neutrinos the expression 

£ M S (mis)- 1 m D + (m s ^- (l + §-) ) , (4.7) 

V m RS V m RS J J 



m u = m D m RS 



which follows from Eq. (1 2 . 1 5 [) after inserting Eq. (14. 5p . The first term on the right-hand 
side represents the well known double see-saw formula. By setting in the suggestive values 
rriD — 100 GeV, Ms — Mpi and ttirs — A^gut — 10 16 GeV, we can generate the correct 
order of magnitude for neutrino masses. The second term in Eq. (14. 7p gives the order 
of magnitude of the very lengthy NLO corrections, which however can be obtained in a 
straightforward manner by inserting Eq. (14.61) in Eq. (I2.15p . For the sake of completeness, 
let us quote the result: 

-2ml = 

m D m ~ R l s M s (rn&g)- 1 (m D m ] D (m^) -1 M* s (mj^) -1 + (last term) r j m R l s M s (ml s )~ X ™d 
+ m T D m^s (m s {m T RS )~ l m D m ] D {m* RS )~ l + (last term) T j {m T RS )~ l m D . 

(4.8) 

The first term of order M| m A D m RS is the leading one for the double see-saw. By setting 
in the suggestive values given above, we realize that the latter gives a contribution of the 
same order, if not larger, than the correction m A D /M R for the type I see-saw formula. It 
is however, with the indicated values of m^, Ms and m^s a negligibly small contribution, 
which may change in other realizations. 

Let us now determine the unitarity violation. From Section [2~731 we know that we can read 
off its amount from the first order expression of B (cf. Eq. ( 12 . 3 1 [) ) . Thus, 

r] ~ — i Mjr, (M^ 1 )* M^ 1 M d i (4.9) 
which has the following explicit form: 



\[m) D {m RS ) 1 M* s {m) RS ) l m RS M s {m T RS ) 1 m D + m j D (m* RS ) 1 {mis) 1 m o) 



[3- (1 + ^i- 

RS \ m RS 



(4.10) 

Again, for the double see-saw the term Msm 2 D /m RS in the second row, giving the order 
of magnitude of rj, is expected to be the dominating one. 

Finally, we mention the possibility of "screening" [22], in which case = em]J s . It 
follows that the zeroth plus first order terms are given as 



m, y = e 2 



Ms- 1 -^ M s [m* s (m^)- 1 + (m^)" 1 M*) M s - e 4 M s . (4.11) 
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One sees that the leading order term has its flavor structure determined by the high 
(possibly Planck) scale physics, while the corrections include additional flavor terms. The 
unitarity violation simplifies to 

V = ~\t 2 (l + M s (m^)" 1 m-i M s ) . (4.12) 

4.2 Inverse See-Saw 

The inverse see-saw [6] is a variant of the double see-saw. The texture of the neutral fermion 
mass matrix is the same as in Eq. ( 14. 2ft . but now obeys the condition Ms <C mo <C wirs- 
In the limit of Ms — > lepton number is conserved and the scenario is natural in the 't 
Hooft sense [23]. It is the preferred scenario to arrange for sizable unitarity violation. A 
recent discussion can be found in Ref. [21], and the results for 77 in this paper agree. 
The calculation of the NLO terms proceeds in the same way as for the double see-saw, 
because we can perform the same replacement as in Eq. H4.3[) : the eigenvalues of M. R 
(which form a Pseudo-Dirac pair) are much larger than the entries in Mj). The effective 
light mass matrix m„ and the unitarity violating parameter 77 look exactly as in Eq. (14. 7p 
and Eq. (]4.10p . respectively. However, the term of order M s m 2 D /m RS is not anymore the 
dominating one, but can be neglected instead. This means that 77 does basically not depend 
on Ms and is given by 

r]---m ] D {m* RS Y 1 {m T RS Y l m D . (4.13) 

With the suggestive values m_o = 100 GeV, m R s — 1 TeV and Ms = 0.1 keV it follows 
that 77 is of order 10~ 2 . The leading term of the NLO correction to the mass matrix is 

ml = -| m T D m RS [m s {m T RS )~ l m D m} D [m\ s y x + (last term) T ) (m^)" 1 m D . (4.14) 

It is of order mj^mj^ Ms and for mo = 100 GeV, m R s — 1 TeV and Ms = 0.1 keV of 
order 10~ 2 eV. In analogy to the case treated in Section I2.2[ the sizable NLO term has 
its origin in the moderate hierarchy of two sub-matrices in the total neutral fermion mass 
matrix. 

4.3 Linear See-Saw 

The linear see-saw mechanism pj arises when the neutral fermion mass matrix has the 
following form: 

/ m T D m T DS \ 
M= m D mis ■ ( 4 - 15 ) 
\ m DS m R s M s J 

Here the non-zero 31 entry mos is assumed to be of weak scale, i.e. of order mo- In the 
following we will assume that niRs is much larger than 777.73 and 771795. Its flavor structure 
may or may not be related to the flavor structure of 777755. The relations between the other 
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block matrices in Ai are that of the double or the inverse see-saw and are given above. We 
can introduce the notation 

M D :=( mD ) and M R := ( ° "ff ) , (4.16) 

V m DS J V mRS Ms J 

and write 

M=( ° M M (417) 

V M D Mr J ^ U > 

The eigenvalues of are much bigger than the entries in Mbj which allows us again to 
apply the method of Section [2] on matrix (I4.17p . The uncorrected low energy mass matrix 
is easily obtained as 

m° = m T D m RS M s (m^ 5 ) _1 m D - [m T D m~ R \ m DS + m T DS {m T RS )~ l m D ] . (4.18) 

Note that if the first term was negligible and m^s oc ttirs, the flavor structure of m u is 
determined by the flavor structure of tub- The unitarity violation rj is again determined 
by Eq. (j2.3ip . There are in total 5 different terms, two of which are the known ones from 
Eq. (I4.10p . and the remaining three are 

V - V douUc ~ \ (m DS (rriRs)' 1 m DS - m ] D (m^) -1 M* (m^)" 1 m^ s m DS 

-™ds (^rs)' 1 m RS M s (mis)' 1 ™d) (4.19) 
= (^(l + ^k + M s m ° + M% , ml , )) ■ 

\ m RS \ m DS ^rriDsmRs m RS m DS ) ) 

As mentioned above, quite often it holds in explicit realizations that rrios — ^ m RSi i n 
which case the first new term in Eq. ( I4.19P is proportional to e 2 1, and the contribution 
to the mass matrix is —e(mo + mj). If we assume that the terms containing Ms are 
absent or sufficiently suppressed, then unitarity violation is determined by terms of order 
( m i>5+ m z))/ m RS- Sizable violation of unitarity could be achieved if mr> or ttids are sizable 
and not much smaller than turs- 



4.4 Singular See-Saw 
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(4.20) 



Cases with a vanishing determinant of Mr are called singular see-saw [TJ, and have re- 
cently received some new attention in the framework of sterile neutrino hints in LSND or 
MiniBooNE data [25]. We will shortly apply our approach to this case now. 
In a three generation framework, the mass matrix is 



M 
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There are two heavy mass states of order M\^, two light states of order m D /Mn and two 
intermediate states of order mp, which form a Pseudo-Dirac pair. Realistic cases with 3 
light neutrinos would require that mo and Mr are 4x4 matrices, the latter having rank 
3. 

We can remove first the heavy states from the discussion by identifying 



/ a x \ 

■ a 2 

• • a 3 

V • • • / 



id 



bi b 2 b 3 

ci c 2 c 3 



Mi 
M 2 



(4.21) 



The low mass states (i.e., the small masses and the Pseudo-Dirac pair) are obtained from 
diagonalizing 
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t 1 ! t>2 

Mi 
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Cl C2 

,2 



^*2 __|_ ^2 

m7 Ml 
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fc 2 h 
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+ M 2 

_|_ C2 C3 
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1,2 .2 
°3 I _£3_ 

Mi "r" M 2 



-a 2 



-a 3 




7 



(4.22) 



The corrections to this matrix can be evaluated using the expression for the NLO term in 
case a triplet is present, see Eq. ( I2.16p . There are two terms, one steming from M.r, the 
other from M^. Their structure is different, the contribution from looks like 



/ * * * \ 

* * * 

* * * 

y * * * o j 



(4.23) 



where the non-zero entries are of order m D /M R . 
ture 

/0 















The contribution from Mi has the struc- 



(4.24) 



where the non-zero entries are of order m D /Mj. 



The relative correction to all entries is 



therefore the same, namely of order m 2 D /M R . As mentioned in Section [2~3| a triplet term 
does not contribute to unitarity violation (as long as mi <C Mr), and t] is determined 
solely by Mb and M.r. The result is 



T] = - 
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N 2 | kil 2 



Mi 



Mi 
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-LS. 
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bx b* 
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leal 2 





~Mj 




Mf 



(4.25) 
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Though its entries are arbitrary, rj is effectively only a 3 x 3 matrix, having no effect for 
the fourth state, which is one of the Pseudo-Dirac states. 

If for instance in the double see-saw of Section 14.11 the matrix rrtRs was singular, we could 
now apply similar steps. After suitable diagonalization of M# in Eq. (14.31) it would (recall 
that M s ^> m RS ) take a form corresponding to diag(0, m 2 RS /M s , m 2 RS /M s , M s , M s , Ms). 
Here the entries are understood as being of order m RS /M s and M s , respectively. Because 
of m 2 RS /Ms ^> vrtY) we can redefine M# as being a diagonal 5x5 matrix of the form 
dia,g(m 2 RS /Ms, m 2 RS /Ms, Ms, Ms, Ms) and follow the procedure of this subsection, finding 
a Pseudo-Dirac pair in the general case etc. Interestingly, if turs was rank 2, we could 
write it as 



The eigenvector of the vanishing eigenvalue of M# (which has rank 5, if Ms is non-singular) 
is proportional to (0, — 1/c, 1, 0, 0, 0) T , i.e. a similar situation as for scaling treated in Section 
13.21 Note that with mus having rank 2, and being rank 5, the full 9x9 mass matrix 
has rank 9; there is no vanishing eigenvalue. We will not discuss the cases of "singular 
double see-saw" or "singular inverse see-saw" any further. 

5 Conclusions and Summary 

With increasing precision in the experimental determination of neutrino mass and lepton 
mixing parameters, care has to be taken in giving theoretical predictions. In the present 
paper we have revisited higher order corrections to the see-saw mechanism. The conven- 
tional type I see-saw, as well as several popular variants were considered, and a strategy to 
determine the next-to-leading order (NLO) terms was developed, based on the well-known 
formalism for the type I see-saw. This can be applied to determine both the NLO terms, as 
well as to obtain the structure of the unitarity violation connected to see-saw mechanisms 
with additional neutral fermions. Table [1] summarizes the structure of the zeroth and 
next-to-leading order terms, as well as of the parameter describing the unitarity violation. 
We have identified situations in which no corrections arise to certain observables, namely 
vanishing neutrino masses or /i-r symmetry. 

While the standard type I see-saw implies insignificant NLO terms, there are cases with 
phenomenologically interesting effects. This occurs when sub-matrices of the complete neu- 
tral fermion mass matrix obey a moderate hierarchy (say, TeV and weak scale). Examples 
are scenarios which explain the smallness of neutrino masses through cancellations of large 
terms, or inverse see-saw frameworks. NLO corrections in the percent regime can arise in 
those cases. 




(4.26) 
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Table 1: See-saw variants and their "typical" orders of magnitude for the zeroth order mass matrix m°, the NLO term 
ml and the unitarity violating parameter rj. 



